FAST TRACK COMMUNICATIONS 



Breaking down the Fermi acceleration with inelastic 
collisions 

O 

O ' Edson D. Leonel 

CN . 

, Departamento de Estatistica, Matematica Aplicada e Computagao - IGCE - 

■ Universidade Estadual Paulista - UNESP - Av. 24A, 1515 ~ Bela Vista - 13.506-900 

- Rio Claro - SP - Brazil 

Abstract. The phenomenon of Fermi acceleration is addressed for a dissipative 
bouncing ball model with external stochastic perturbation. It is shown that the 
introduction of energy dissipation (inelastic collisions of the particle with the moving 
wall) is a sufficient condition to break down the process of Fermi acceleration. The 
phase transition from bounded to unbounded energy growth in the limit of vanishing 
dissipation is characterized. 
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The phenomenon of Fermi acceleration (FA) is a process in which a classical particle 
acquires unbounded energy from collisions with a heavy moving wall [1]. This 
0> . phenomenon was originally proposed by Enrico Fermi [2] as a possible explanation for the 
origin of the large energies of cosmic particles. His original model was later modified by 
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several investigators and applied in several fields of physics, including plasma physics [3], 
astrophysics [H |5] , atomic physics [6] , optics [3 El [9] and the well known time-dependent 
billiard problems [ini[ll]. Since the seminal paper of Hammersley [12], it is known that 
the particle's average energy grows with time when there is a random perturbation at 
each impact with the wall. This result was also confirmed for a stochastic version of the 
one-dimensional bouncing ball model (a classical particle confined in and hitting two 
rigid walls; one of them with a fixed position and the other one with periodic movement) 
under the framework of random shift perturbation [131 [13! • 

One of the most important questions on Fermi acceleration is whether it can result 
from the nonlinear dynamics in the absence of a random component . The answer to this 
question depends on the model under consideration. For example, for a bouncer model (a 
particle hitting a periodically moving platform in the presence of a constant gravitational 
field) , there are specific ranges of control parameters and initial conditions that lead to 
Fermi acceleration [15]; FA occurs when there are no invariant spanning curves on the 
phase space limiting the chaotic sea and, consequently, the particle's energy gain [T3lfTi]. 
For two-dimensional time-dependent billiards (billiards with moving boundaries), the 



Breaking down the Fermi acceleration with inelastic collisions 



2 



emergence of FA depends on the type of phase space of the corresponding static version 
of the problem [16]. Therefore, as conjectured by Loskutov and collaborators [16j, the 
chaotic dynamics of a particle for static boundary is a sufficient condition to produce 
FA if a time perturbation in the boundary is introduced. 

A second and very important question is: in the classical billiard problems where 
FA is present, are inelastic collisions of the particle with the boundaries a sufficient 
condition to suppress the unlimited energy gain? 

In this Letter, we consider the one- dimensional Fermi accelerator model under 
stochastic perturbation and we seek to understand and describe a mechanism to suppress 
the FA. The model consists of a classical particle confined to bounce between two rigid 
and infinitely heavy walls. One of them is fixed while the other one moves randomly 
with dimensionless amplitude of motion e \T7]. It is assumed that collisions with the 
moving wall are inelastic, so that the particle experiences a fractional loss of energy 
upon each collision. A restitution coefficient a G [0, 1] controls the strength of the 
dissipation. For a = 1, all collisions are elastic and therefore FA is observed [I5]. On 
the other hand, for a < 1 the model is dissipative and, as shown here, there is a bound 
to the energy growth. In other words, inelastic collisions break down the FA process. 
A phase transition from bounded to unbounded energy growth is observed when the 
control parameter a approaches the unity. Here we investigate this phase transition 
by studying the limit a ^ 1. The present approach can be useful as a mechanism to 
subdue the phenomenon of FA in time- dependent billiard problems. 

The dynamics of the problem is given by a two-dimensional, nonlinear mapping for 
the particle velocity V and the time t at each impact of the particle with the moving wall. 
We investigate in this Letter a simplified version of the model [HI [HI [THl [19], [20] which 
speeds up the numerical simulations significantly without affecting the universality class. 
However, similar results would indeed be obtained for the full model. The simplified 
version assumes that both walls are fixed but that, when the particle hits one of 
them, it exchanges energy and momentum as if the wall were moving randomly. Thus, 
considering dimensionless variables and taking into account the inelastic collisions, the 
mapping that describes the dynamics of the model is 



where n is the iteration number and Z{n) G [0, 27r) corresponds to a random shift in the 
phase of the moving wall. It is clear that the model has two control parameters, namely 
e and a, whose effect must be considered. 

The most natural observable in problems involving FA is the average velocity, which 
is calculated here in two steps. The first step consists in averaging the velocity over the 
orbit for a single initial condition. It is defined as 
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Figure 1. (Color online) (a) Behavior of V as a function of the iteration number n 
for different values of the control parameters e and a. (b) Curves shown in (a) after 
the change of variables n — > ne^. 



where the index j refers to the jth iteration of the sample i. The second step is to 
calculate the average over an ensemble of M different initial conditions 

1 M 

The behavior of the average velocity V for different values of the control parameters 
is shown in Fig. [H Note that the growth of V with n is described by different curves 
for different values of e (Fig. [It^a)). However, a transformation n —>■ ne^ coalesces all 
curves, so that they grow together for small n, as seen in Fig. [U^b). Also, instead of 
studying the behavior of as a function of the damping coefficient a, we adopt the 
variable (1 — a) in order to bring the transition to the origin. Such a transformation 
improves visualization in log-log plots. 

Let us now discuss the behavior observed in Fig. [1] It is easy to see that all curves 
start growing together for small iteration numbers and then they bend towards a regime 
of convergency. Such a regime is marked by a constant plateau for the average velocity. 
The change from growth to saturation is characterized by a typical crossover iteration 
number n^. Based on the results seen in Fig. [H one concludes that the velocity grows 
as a power law of the type V oc {ne^Y for n <^ n^- For large iteration numbers, the 
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saturation velocity for a fixed damping coefficient a is given by V^sat e"^. On the 
other hand, for a fixed e, the values obtained are Vsat oc (1 — a)"^. The typical crossover 
iteration number that marks the transition from growth to the saturation is assumed 
to be of the following form: (i) for a constant a, n^e^ oc e^^ and; (ii) for a fixed e, 
Hxe^ oc (1 — aY"^. These initial assumptions allow us to propose the following scaling 
hypotheses: 

• For small iteration numbers (n <^ n^,), the average velocity grows as 

V oc [neY , (4) 
where /? is a critical exponent. 

• For large iteration numbers (n ^ n^;), the constant plateau of the velocity is given 
by 

FsatOce-ni-«r , (5) 
where both ai and ^2 are critical exponents. 

• The crossover iteration number is given by 

n^e^ o^e'\l-ay^ , (6) 
with zi and Z2 being the dynamical exponents. 

These three hypotheses allow us to formally describe the average velocity using a 
scaling function of the type 

F(r^e^ e, (1 - a)) = lV{tne^, fe, r{l - a)) , (7) 

where / is a scaling factor and a, h and c are the scaling exponents which must be related 
to the critical exponents /9, ai and ^2, Zi and Z2. Since I is a scaling factor, we can 
choose / = [ne^]~^/°. Using this expression for /, Eq. ([71) is rewritten as 

F(ne^ e, (1 - a)) = [neX-^V,{[ne^]--^e, [ne^Y^l - a)) , (8) 

where the function Vi is assumed to be constant for n ^ n^. Comparing equations ([8]) 
and (jlj), we obtain (3 = — 1/a. After conducting extensive numerical simulations, it was 
found that (3 = 0.494(1) = 0.5, which implies a = —2. Let us now consider the case 
n ^ Ux- This case allows us to choose two distinct values for /, namely: (a) / = t"^^'^ 
and (b) / = (1 - a)"^/'. For ), the scaling function is rewritten as 

F(ne2,e, (1- a)) = e-^y2(e"^^e',e"^ (!-«)) , (9) 

with the function V2 being constant for n ^ and a constant. A comparison of 
equations iQ and ([5]) furnishes ai = —1/6. After fitting a power law to data on the 
plot Fsat X e, we find ai = 1.0011(2) = 1 (see Fig. [2](a)), which provides b = -1. We 
now have to consider case (b), i.e., / = (1 — a)^^^'^. Using this scaling factor, Eq. (J?!) is 
given by 

F(ne^ e, (1 - a)) = (1 - a)---V3{{l - a)"^ne^ (1 - a)---e) (10) 
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where we assume that V3 is constant for n ^ and e constant. Comparing equations 
ffTOj) and ([5]), it is easy to see that 02 = —1/c. A power-law fitting to the data for 
Fsat X (1 - a) gives that 0^2 = -0.5002(1) = -0.5, yielding c = 2. 

Considering the different expressions for the scaling factor / obtained for n <^ 
and case (a) of n n^., we obtain 

ne^ = et . (11) 

A comparison of equations ffTTl) and ([6]) provides zi = a/b. Fitting a power law to the 
data on the plot n^e^ x e furnishes Zi = 2.001(2), as shown in Fig. El^a). This result is in 
agreement with the ratio a/b = 2. As a next step, we consider the case of the different 
expressions for the scaling factor / obtained for n <^ and case (b) of n ^ n^. Such 
procedure gives 

ne^ = (1 -a)t . (12) 

Comparing equations ffT^ and we obtain that Z2 = a/c. A power-law fitting to the 
data on the plot n^e^ x (1 - a) (Fig. Mh)), with e = 10"^ gives Z2 = -1.018(2). This 
result is also in good agreement with the result obtained from a/c = —1. 
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Figure 2. (Color online) (a) Plot ofVsat x e. A power law fitting gives ai = 1.0011(2). 
(b) Plot of V^sat X (1 ^ ct)- The exponent obtained via a power-law fitting is 
a2 = -0.5002(1). 
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Let us now discuss the consequences of the critical exponents 02 and Z2 on the FA. 
It is clear that, in the limit a 1, both Eqs. ([5]) and ([6]) diverge. Therefore, in the 
limit of vanishing dissipation, the average velocity V and the crossover iteration number 
rix diverge. This is a clear signature of a phase transition from bounded to unbounded 
energy growth. Finally, in order to check the validity of the scaling hypotheses, we can 
now proceed to collapse all the curves onto a single and universal plot, as shown in Fig. 

SI 

In summary, the problem of a stochastic bouncing ball model with inelastic 
collisions is addressed. It is shown that, when dissipation is present, the average 
velocity of a classical particle grows with time and then reaches a regime of saturation. 
However, in the limit of vanishing dissipation, the phenomenon of Fermi acceleration 
is recovered. Similar results were also observed for a dissipative and deterministic 
bouncer model [211 ES]. These results allow us to propose the following conjecture: 
For one- dimensional billiard problems that show unlimited energy growth for both their 
deterministic and stochastic dynamics, the introduction of inelastic collision in the 
boundaries is a sufficient condition to break down the phenomenon of Fermi acceleration. 
Such phenomenon is also expected to be observed in two-dimensional, time-varying 
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Figure 3. (Color online) (a) Plot of Uxe^ x e. After fitting a power law we obtain that 
zi — 2.001(2). (b) Plot of n^e'^ x (1 — a) for a fixed e = 10"^. Fitting a power law, we 
find Z2 = -1.018(2). 
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Figure 4. (Color online) Collapse of the curves shown in Fig. [TJb) onto a single and 
universal plot. 



billiard problems since the FA mechanism is the same as that of the one-dimensional 
case. 
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